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ABSTRACT 

An  optimization  approach  to  prescribe  open  boundary  conditions  (OBCs)  is  pro¬ 
posed.  Open  boundary  conditions  are  chosen  by  providing  the  "best  fit  to 
available  observations  on  the  open  boundary  and  to  the  energy  flux  through  the 
open  boundary.  It  is  shown  that  many  "local”  open  boundary  conditions  widely 
used  in  the  oceanographic  community  are  special  cases  of  the  derived  optimized 
OBCs.  Numerically,  the  proposed  methods  consist  of  the  integration  of  the  gov¬ 
erning  equations  and  solving  the  optimization  problem  only  on  the  open  bound¬ 
ary  for  each  time  step.  Results  of  tidal  simulations  for  a  channel  and  the 
Louisiana-Texas  shelf  are  presented  and  discussed.  Schemes  for  coupling  limited  area 
coastal  models  with  basin  scale  coarse  models  by  using  optimized  OBCs  are  pro¬ 
posed.  Preliminary  simulation’  of  coupling  a  fine  resolution  model  of  the  north¬ 
ern  part  of  the  Adriatic  Sea  and  a  coarse  resolution  model  of  the  Mediterranean 
Sea  are  presented  and  discussed.  The  use  of  optimized  OBCs  is  shown  to  be 
better  in  predicting  tidal  amplitudes  and  phases  compared  to  the  use  of  "local” 
OBCs. 


1.  INTRODUCTION 

For  operational  needs  of  the  Navy,  coastal 
nowcast/forecast  models  will  include  large  open 
boundaries  in  the  domain  of  interest  because 
there  is  no  single  numerical  ocean  model  which 
can  provide  the  necessary  continuous  resolution 
from  the  open  ocean  to  the  coastal  area.  The 
treatment  of  open  boundaries  is  one  of  the  most 
interesting  problems  to  be  solved  while  mod¬ 
eling  oceanic  phenomenon,  especially  in  finite 
ocean  coastal  areas.  In  most  ocean  models,  open 
boundary  conditions  (OBCs)  are  chosen  locally, 
i.e.,  depending  on  the  solution  of  the  governing 
equations  near  the  boundary.  Many  approaches 
of  the  local  type  have  been  developed  (Reid  and 
Bodine,  1968;  Orlanski,  1976;  Chapman,  1985; 


Blumberg  and  Kantha,  1985;  Oey  and  Chen, 
1992  (originally  due  to  Flather,  1976),  etc.).  The 
results  of  numerical  studies  show  that  the  applica¬ 
tion  of  many  local-type  boundary  conditions  re¬ 
produce  the  modeled  physical  phenomenon  and 
work  for  most  practical  purposes.  However,  it 
is  known  (Bennett,  1992;  Oliger  and  Sundstrom, 
1978)  that  the  local  treatment  of  open  boundaries 
for  primitive  equation  models  is  an  ill-posed  prob¬ 
lem  in  that  it  is  difficult  to  prove  that  a  unique 
solution  exists  that  is  continuously  dependent  on 
available  observations.  Some  researchers  use  the 
inverse  approach  to  the  modeling  of  open  bound¬ 
ary  conditions.  The  open  boundary  conditions 
are  chosen  in  such  a  way  as  to  simultaneously 
provide  the  "best”  fit  to  the  governing  equations 
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and  observations.  The  "best"  fit  means  the  min¬ 
imization  of  the  norm  of  the  deviation  between 
model  results  and  observations.  Thus,  the  inte¬ 
rior  solution  and  available  observations  are  used 
to  calculate  variables  on  the  open  boundary.  This 
approach  has  been  applied  to  atmospheric  and 
ocean  circulation  problems  (e.g.,  Bennett,  1992; 
Zou  et  al.,  1993;  Bogden  et  al.,  1994;  Seiler, 
1993).  The  most  popular  algorithm  for  solving  an 
inverse  problem  is  an  adjoint  method  in  which  the 
initial  problem  of  circulation  with  open  boundary 
conditions  is  reduced  to  integrating  the  governing 
equations  and  an  equation  for  the  adjoint  vari¬ 
able  forward  and  backward  in  time  (Zou  et  al., 
1993).  Although  the  inverse  approach  leads  to 
weil-posedness,  it  suffers  from  a  few  drawbacks 
that  may  restrict  its  use:  requirements  for  large 
amounts  of  computer  time  and  memory  and  the 
problem  of  the  stable  integration  of  the  adjoint 
equation.  An  even  more  serious  problem  in  using 
an  inverse  method  is  not  knowing,  a  priori,  the 
physical  properties  of  the  solution.  In  spite  of  its 
“best'' ,  the  solution  may  be  physically  wrong. 

There  are  two  major  aspects  in  developing 
OBCs  for  coastal  models.  The  first  one  is  to  cor¬ 
rectly  formulate  OBCs  mathematically  and  phys¬ 
ically  for  the  purpose  of  reproducing  the  modeled 
phenomenon.  The  second  one  is  to  develop  nec¬ 
essary  techniques  for  OBCs  initialization  by  the 
required  data.  These  data  can  be  derived  from  a) 
available  observations,  b)  coarse  resolution,  basin 
scale  model  simulations,  or  c)  a  combination  of 
a)  and  b).  Initialization  can  be  accomplished  by 
the  development  of  the  technology  for  coupling  a 
coastal  model  and  a  basin  scale  model  with  the 
use  of  data  assimilation  on  the  open  boundary. 

In  this  paper,  we  propose  methods  for  mod¬ 
eling  open  boundary  conditions  which  are  based 
on  the  integration  of  the  governing  equations  for¬ 
ward  in  time  and  calculating  values  of  variables 
on  the  open  boundary  via  a  specific  inverse  prob¬ 
lem  that  provides  the  "best”  fit  to  available  ob¬ 
servations  on  the  open  boundary  and  to  the  en¬ 
ergy  flux  through  the  open  boundary  (Shulman 
and  Lewis,  1994).  In  this  way,  we  avoid  the  lo¬ 
cal  treatment  of  the  open  boundary  conditions. 
All  observations  on  the  open  boundary  and  its  in¬ 
terior  vicinity  are  used  in  determining  the  open 
boundary  conditions  for  each  particular  point  on 
the  boundary.  Numerically,  the  proposed  meth¬ 
ods  consist  of  the  integration  of  the  governing 
equations  and  solving  tne  optimization  problem 
for  each  time  step.  We  call  these  open  bound¬ 


ary  conditions  optimized  OBCs.  We  show  that 
some  of  the  well  known  “local”  long-wave  radi¬ 
ation  boundary  conditions  which  are  commonly 
used  in  ocean  modeling  are  special  cases  of  the 
optimized  OBCs  derived  using  this  approach.  The 
derivations  presented  here  suggest  methods  for 
the  generation  of  new  boundary  conditions  based 
on  tne  requirements  of  the  modeled  phenomena. 

Most  of  the  approaches  to  the  coupling  of 
basin  scale  and  coastal  models  can  be  divided 
into  two  groups  (Spall  and  Holland,  1991;  Oey 
and  Chen,  1992):  "one-way”  and  “two-way”  cou¬ 
pling.  In  "one-way”  coupling,  boundary  data  for 
a  limited  area  model  (LAM)  are  prescribed  by  us¬ 
ing  data  from  a  basin  model  (usually  a  coarser 
resolution,  larger  domain  model).  In  “two-way” 
coupling,  the  evolution  within  the  LAM  influences 
the  evolution  within  the  coarser  model.  We  dis¬ 
cuss  an  investigation  of  "one-way”  coupling  with 
the  use  of  our  techniques  for  the  specification  of 
OBCs. 


2.  DERIVATION  OF  THE  OPTIMIZATION 
PROBLEM 

Let  Pt  be  the  estimated  energy  flux  through 
an  open  boundary,  Mt  be  the  estimated  mass  flux 
through  the  same  open  boundary,  Ft  be  the  esti¬ 
mated  momentum  flux  through  the  open  bound¬ 
ary,  and  vector  X  be  the  variables  which  we 
should  specify  on  the  open  boundary  (sea  sur¬ 
face  height,  velocity,  temperature,  salinity,  etc.). 
Pt,  Mt,  Ft  can  be  estimated  by  using  model  val¬ 
ues  from  the  interior  solution  or  from  simulations 
from  another  model  (for  example,  a  larger  do¬ 
main,  coarse  resolution  model).  Let  J(X)  be 
some  function  of  X,  J(X )  >  0.  We  can  include 
available  observations  ana  a  priori  information 
in  the  specification  of  J(X).  Basically,  function 
J(X)  represents  the  difference  between  the  val¬ 
ues  of  model  variables  and  observations  on  the 
open  boundary.  For  each  time  step  we  choose 
open  boundary  conditions  X  from  the  following 
inverse  problem: 


min  J(  X), 

(1) 

A(X)  =  Pt , 

(2) 

B{X)  =  Mu 

(3) 

C(X)  =  Ft, 

(4) 

> 


» 


I 
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where  A,B,C  are  operators  for  calculating  en¬ 
ergy,  mass  and  momentum  fluxes.  Thus,  we 
choose  open  boundary  conditions  by  provid¬ 
ing  the  "best"  fit  to  available  observations 
on  the  open  boundary  and  to  the  estimated 
energy,  mass,  and  momentum  fluxes  through 
the  open  boundary.  Solving  (1W4)  is  a 
very  complicated  problem,  and  different  ap¬ 
proaches  can  be  used.  They  all  have  to  take 
into  account  that  the  estimated  Pt,  Mt,  Ft, 
observations,  and  a  priori  information  in  J(X) 
may  contain  errors. 

In  this  paper  we  only  consider  the  solution 
of  (l)-(2)  for  the  vertically-averaged,  hydrostatic 
equations.  We  introduce  the  notation 


Pt  = -g  J^HrjUnds.  (5) 

where  S  is  the  open  boundary  of  the  model  do¬ 
main  D,  u„  is  the  vertically-averaged  outward 
normal  velocity,  rj  is  the  sea  surface  deviation, 
H  is  the  depth,  and  g  is  the  gravitational  con¬ 
stant.  The  term  Pt  can  be  interpreted  as  the  flux 
of  total  energy  penetrating  the  open  boundary 
5  for  a  model  based  on  the  shallow  water  equa¬ 
tions,  or  as  a  part  of  the  pressure  force  work  for  a 
vertically-averaged,  hydrostatic  model.  Suppose 
we  estimated  the  energy  flux  Pt.  We  can  intro¬ 
duce  a  criteria  to  choose  one  specification  of  un 
and  77  which  satisfies  (5).  Consider  the  following 
optimization  problem: 

min  J(rj),  (6) 


constrained  by  the  flux  of  energy  through  the 
open  boundary: 

Pt  =  ~g  J  Hr)unds.  (7) 

Here  J(rj)  is  an  objective  functional  depending 
on  77,  such  that  J{r 7)  >  0.  Using  the  Lagrangian 
method  (Fletcher,  1987)  to  solve  the  problem  (6)- 
(7),  we  minimize: 

min  [J(r})  —  \t(Pt  +  g  /  Hrjunds )] 

Js 


where  \t  is  a  constant  (the  Lagrangian  multi¬ 
plier).  Thus,  the  solution  of  the  problem  satisfies 
the  following  condition: 


dJ 

- - A  tgHun  =  0. 

dr] 


(8) 


it  is  important  to  note  that  the  constant  A,  mea¬ 
sures  the  rate  of  change  in  the  function  J(g)  due 
to  changes  in  Pt.  Let  et  be  the  perturbation  of  Pt  : 
P*t  =  Pt  +  €<.  It  can  be  shown  that  (Fletcher, 
1987): 


At 


dJ_ 

de,' 


(9) 


which  provides  us  with  a  direct  relationship  be¬ 
tween  the  various  terms  of  the  problem. 


3.  SOME  EXAMPLES  OF  OPTIMIZED  OBCs 

In  the  following,  we  will  show  that  several  of 
the  open  boundary  conditions  used  in  numerical 
ocean  modeling  can  be  related  to  the  optimiza¬ 
tion  approach  of  the  last  section.  In  each  case,  a 
specific  objective  functional,  J,  is  minimized. 

3.1  Simple  Longwave  Radiation 

Let  us  consider  the  following  problem, 

min(7i  =  |  [  y/gHg2  ds),  (10) 

1  *  Js 

with  condition  (7).  In  this  case,  (8)  becomes 


Note  that  (11)  is  the  commonly-used  long-wave 
radiation  condition  with  a  "tuning"  coefficient  At 
which  is  time  dependent.  With  this  coefficient, 
the  relationship  given  by  (11)  minimizes  the  func¬ 
tional  J(rj)  given  in  (10)  (in  essence,  the  square 
of  the  sea  surface  height  on  the  open  boundary). 
If  At  =  1  (the  common  radiation  condition),  then 
according  to  (9)  we  have  dJ/det  =  1.  Under  such 
conditions,  we  are  specifying  that  the  change  in 
the  flux  of  energy  penetrating  the  open  bound¬ 
ary  will  be  equal  to  the  change  of  the  functional 
specified  in  (10). 

3.2  The  Reid  and  Bodine  Boundary 
Condition 

There  are  situations  when  some  extraneous 
information  on  the  open  boundary  in  the  form  of 
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sea  surface  height,  t)t,  is  available.  This  infor¬ 
mation  could  be  direct  observations  or  output  of 
another  numerical  model.  In  this  case,  we  can 
specify  the  objective  functional  to  be: 

min(J2  =  |  /  VgH(v  ~  Vt)2  ds),  (12) 
i  *  Js 

where  the  model  height  is  fitted  to  the  known 
function  in  the  least  squares  sense.  In  this  case, 
(8)  becomes 


Condition  (16)  is  a  modification  of  the  boundary 
condition  employed  by  Oey  and  Chen  (1992)  with 
r]*T  calculated  from  (15).  The  boundary  condi¬ 
tions  given  by  (13)  and  (15)-(16)  both  minimize 
the  difference  between  the  modeled  and  observed 
surface  heights  [equation  (12)].  The  difference  is 
that  (15W16)  allows  for  tne  inclusion  of  currents 
observer/  at  tne  boundary. 

3.3  The  Blumberg  and  Kantha  Boundary 
Condition 

We  specify  the  objective  functional  to  be: 


(v-it)  [7 

u"  -  — —y  h  (13) 

Note  that  condition  (13)  is  a  modified  version 
of  that  used  by  Reid  and  Bodine  (1968),  which  is 
often  used  to  force  regional  tidal  models  while  still 
allowing  for  the  radiation  of  longwave  energy  from 
the  interior  domain  of  the  model  (Lewis  et  al . , 
1992).  In  this  case,  the  relationship  given  by  (13) 
minimizes  the  deviation  of  the  sea  surface  height 
from  the  known  function  rjr  (e.g.,  observed  sea 
surface  height  variations).  When  At  =  1,  the 
change  in  the  flux  of  the  energy  penetrating  the 
boundary  S  is  equal  to  the  change  of  the  deviation 
t]  from  tjt- 

Oey  and  Chen  (1992)  incorporated  current 
velocity  information,  UT(s,t\  into  their  open 
boundary  condition  (originally  due  to  Flather, 
1976).  To  do  so  in  our  formulation,  we  rewrite 
(13)  in  the  following  form: 


min  J3  = 

5  Is  +  ^~rfT^dt) 

where  T  is  a  known  constant.  The  minimization 
of  J3  subject  to  the  energy  flux  constraint  leads 
to: 


T1-Xtj^gUn  =  ~fj0  (r,~TlT)dt‘  (18) 


By  differentiation  of  the  equation  (18)  with  re¬ 
spect  to  t,  and  after  simple  transformations  we 
can  get  the  following  (Shulman  and  Lewis,  1994): 


d(\t) 

dt 


-  (j)T  +  =  At(un  -  U„,t), 


(14) 

where  un>T  is  the  outer  normal  component  of  the 
known  velocity  ut(s,  t).  Let  us  introduce  the  fol¬ 
lowing  notation: 


V*T  —  VT  +  At 


(15) 


- ={t]-t)T )•  (19) 

The  condition  specified  in  (19)  is  a  variation  of 
the  boundary  condition  developed  by  Blumberg 
and  Kantha  (1985)  (if  At  =  1  for  all  t).  According 
to  (17),  condition  (19)  minimizes  the  square  cf  77 
ana  the  deviation  of  77  from  777-  over  time  along 
the  open  boundary.  If  T  =  t2 ,  we  are  minimizing 
the  square  of  77  and  the  average  deviation  of  77 
from  tjt  over  time. 


By  substituting  (15)  into  (14),  we  can  get  the 
following: 


3.4  Some  Enhancements 

A  few  simple  enhancements  of  the  objective 
functionals  can  be  considered  to  derive  different 
open  boundary  conditions.  Consider 

min  J4  =  0\J\  +  02  J2  = 
v 
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|  /  VgBlPrf  +  Hi  -  1Tf]ds,  (20) 

again  with  constraint  (7).  For  this  case,  (8)  be¬ 
comes 


We  assume  that  Pt  and  un<t  contain  some  er¬ 
rors  in  the  estimates  of  the  energy  flux  and  the 
normal  velocity.  We  employ  the  regularization 
method  (Sabatier,  1987;  Parker,  1994)  for  solv¬ 
ing  the  problem: 


@\T)  -  A t 


-&(»7  -  Vt)- 


If  =  1  and  fl 2  =  we  have: 


1  f 

minl ^{Pt+9  Hriun  t  ds)  +  7  J],  (24) 

^2  Js 

where  7  is  a  parameter  of  regularization.  The 
solution  of  (24)  gives  the  following  expression  for 
\t  in  (13): 


n 


-=(rj -tit)- 


(21) 


Pt  +9  Js  Hr)t  TUn,t  ds 
gifsHiun,t2ds  +  7' 


(25) 


This  formulation  provides  for  the  direct  specifica¬ 
tion  of  irn  at  the  boundary. 


4.  DISCUSSION 

The  optimized  open  boundary  conditions 
(11),  (IS),  (16),  (19),  and  (21)  are  not  local 
point  wise  OBCs,  because  for  calculating  At  we 
use  data  all  along  the  open  boundary  and  its  in¬ 
terior  vicinity.  There  are  many  possible  numerical 
approaches  to  implementing  the  proposed  bound¬ 
ary  conditions.  These  approaches  depend  on  the 
numerics  of  the  hydrodynamical  model.  Here  we 
discuss  one  of  these  approaches  for  the  optimized 
Reid  and  Bodine  boundary  condition  (13)  imple¬ 
mented  for  c  vertically-averaged  model.  Suppose 
rjt  and  un  t  are  the  sea  surface  elevation  and  ve¬ 
locity  at  time  t,  and  a  model  uses  two  previous 
time  steps  for  calculating  variables  for  the  t+1 
time  step.  As  always,  the  sea  surface  elevations 
77<+i  are  calculated  from  the  continuity  equation, 
and  the  velocities  itt+1  are  calculated  from  the 
momentum  equation.  On  the  open  boundary, 
the  following  numerical  scheme  was  implemented 
based  on  our  optimization  approach.  By  using 
values  of  r)t  and  un>t,  we  can  find  an  estimate  of 
Pt  from  (5).  We  have  the  following  optimization 
problem  for  time  t: 

min(J=  ^  /  y/gH{q  -  riTj)2  ds),  (22) 
Pt  = -g  J  Hrjun,tds.  (23) 


After  the  calculation  of  \t,  we  rearrange  (13)  to 
solve  for  a  corrected  value  of  sea  surface  elevation 
at  time  t  using  the  value  of  un  calculated  for 
time  t.  This  corrected  value  of  the  sea  surface 
elevation  at  time  t  is  used  for  calculating  the  value 
of  «n,t+ 1 

u„,t+i  =  (r]t  -  rlT,i)\j^-  (26) 

Because  we  do  not  know  the  norms  of  the  errors 
in  the  estimates  of  Pt  and  unj(,  the  parameter  7 
is  chosen  in  such  a  way  that  the  two  terms  in  (24) 
are  equal. 

5.  NUMERICAL  EXPERIMENTS 

We  performed  simulations  for  a  flat  bottom, 
frictionless  channel  which  is  closed  at  one  end. 
The  channel  was  forced  at  the  other  end  by  sur¬ 
face  oscillations  at  the  M2  tidal  frequency  with 
an  amplitude  of  1  m.  The  length  of  the  channel 
was  355  km.  There  are  grid  points  along  the 
channel  (the  23rd  is  a  wall).  We  performed  simu¬ 
lations  for  three  different  values  of  the  channel’s 
depth:  255  m,  1000  m,  and  50  m.  Figs.  1-3  show 
the  amplitudes  and  phases  of  the  standing  waves 
based  on  the  analytical  solution,  a  model  simula¬ 
tion  using  the  Reid  and  Bodine  boundary  condi¬ 
tion,  and  a  model  simulation  using  the  optimized 
Reid  and  Bodine  formulation.  The  use  of  the  op¬ 
timized  Reid  and  Bodine  formulation  is  seen  to 
eliminate  almost  entirely  errors  in  the  predicted 
amplitudes  while  more  than  halving  the  errors  tn 
the  predicted  phases. 
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We  performed  simulations  for  the  M2  tide 
over  the  Louisiana-Texas  shelf  (LATEX)  area 
(Fig.  4).  The  model  domain  has  southern  and 
eastern  open  boundaries,  and  forcing  at  these 
boundaries  was  determined  from  the  gobal  tide 
model  of  Schwiderski  (1983).  Tables  1  and  2 
contain  the  results  of  the  simulations  performed 
with  both  the  Reid  and  Bodine  and  optimized 
Reid  and  Bodine  conditions  [calculating  At  from 
(25)].  The  use  of  the  optimization  approach  gives 
a  better  prediction  of  the  amplitudes  and  phases 
at  Galveston,  Texas;  Port  Aransas,  Texas;  and  the 
weak  tide  at  South  West  Pass,  Louisiana.  Fig.  5 
shows  a  comparison  between  the  tidal  rorcing  at 
the  open  boundaries  and  the  tides  predicted  at 
the  model  grid  cells  next  to  the  open  boundaries. 
These  results  indicate  that  the  modified  Reid  and 
Bodine  formulation  provides  a  good  estimate  of 
the  tidal  amplitudes,  with  the  largest  deviations 
occurring  on  the  eastern  open  boundary  along  the 
shelf  slope.  However,  the  minimization  process 
for  this  simulation  shows  a  clear  bias  in  overesti¬ 
mating  the  phases,  with  the  largest  overestimates 
occurring  along  the  southern  open  boundary. 

We  also  conducted  an  experiment  where 
the  southern  and  eastern  open  boundaries  were 
treated  separately  in  our  optimization  approach. 
Different  values  of  A<  were  calculated  for  each 
open  boundary.  Thus,  the  technique  would  min¬ 
imize  differences  along  the  eastern  open  bound¬ 
ary  independently  of  those  differences  along  the 
longer  southern  open  boundary.  The  results  of 
the  simulation  are  presented  in  Table  3.  In  this 
case,  the  amplitudes  were  overestimated  and  the 
phases  were  underestimated.  As  we  mentioned 
in  the  introduction,  the  original  Reid  and  Bod¬ 
ine  condition  is  a  local  treatment  approach  to 
the  specification  of  the  open  boundary  condi¬ 
tion:  only  neighborhood  data  are  used  in  deter¬ 
mining  the  boundary  condition  for  a  grid  cell  on 
the  boundary.  In  our  approach,  all  data  along 
and  near  the  boundary  are  used  in  calculating  the 
boundary  values  for  each  grid  point  on  the  bound¬ 
ary,  and  this  approach  results  in  a  substantial  im¬ 
provement  in  predicting  amplitudes  and  phases. 
However,  the  separate  treatment  of  the  south¬ 
ern  and  eastern  boundaries  gives  inferior  results, 
likely  due  to  its  more  “local"  treatment  of  the 
boundaries. 


6.  COUPLING  TECHNOLOGY  AND 
NUMERICAL  EXPERIMENTS 

Conditions  (13),  (16),  (19)  and  (21)  can  be 
considered  methods  of  boundary  value  relaxation 
toward  either  observational  data  or  results  from  a 
coarser,  basin  scaie  model  simulation.  These  con¬ 
ditions  have  coefficients  of  relaxation  that  change 
in  time  and  provide  the  adaptation  of  the  bound¬ 
ary  conditions  to  the  change  in  the  energy  flux 
through  the  open  boundary.  We  consider  two 
different  schemes  for  coupling  basin  and  coastal 
models  with  the  use  of  data  assimilation  on  the 
open  boundary: 

1.  providing  the  "best"  fit  to  the  data  on  the 
open  boundary  obtained  from  observations 
and/or  larger  domain,  basin  model  output 
and  to  the  energy  flux  generated  by  an  inte¬ 
rior  solution  of  a  limited  area  model  (LAM), 

2.  providing  the  "best"  fit  to  the  data  on  the 
open  boundary  obtained  from  observations 
and/or  larger  domain,  basin  model  output 
and  to  the  energy  flux  calculated  from  the 
basin  model. 

We  tested  the  first  method  of  coupling  by 
conducting  M2  tidal  simulations  in  the  Mediter¬ 
ranean  Sea  using  the  explicitly  formulated,  sigma 
coordinate  version  of  the  Princeton  ocean  model 
(Blumberg  and  Mellor,  1987).  Two  orthogonal 
curvilinear  grids  were  used  which  cover  the  entire 
Mediterranean  Sea.  The  first  grid  had  a  relatively 
fine  mesh,  with  a  grid  size  ranging  from  4.3  km  to 
12.3  km  and  441x141  grid  points  in  the  horizon¬ 
tal.  The  second  grid  had  a  relatively  coarse  mesh, 
with  a  grid  size  ranging  from  3.5  km  to  62  km  hor¬ 
izontal.  All  simulations  were  performed  with  tne 
bathymetry  based  on  a  2”x2”  data  base  devel¬ 
oped  by  the  Naval  Oceanographic  Office  (NAV- 
OCEANO).  Direct  gravitational  M2  forcing,  tidal 
forcing  at  the  western  open  boundary,  and  details 
about  boundary  conditions  and  bottom  friction 
coefficients  can  be  found  in  Lewis  et  al.  (1994). 

The  results  of  the  M2  tidal  simulations  in 
the  region  of  the  northern  Adriatic  Sea  using  the 
finer  grid  model  are  shown  in  Fig.  6.  One  can 
see  a  good  agreement  between  the  model  and  the 
observational  data  from  six  stations:  Porto  Pi- 
ave  Vecchia,  45.29°N,  12.34°E;  Rovinj,  45.05°N, 
13.38°E;  Porto  Corsini,  44.30°N,  12.17°E;  Pesaro 
43.55°N,  12.55°E;  Ancona,  43.37°N,  13.30°E; 
and  Pula,  44.52°N,  13.50°E.  To  quantify 
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Fig.  5.  Amplitude  (left)  and  phase  (right)  of  the  M2  tide  along  the  open  boundaries  of  the  LATEX 
shelf  model  using  the  modified  Reid  and  Bodine  open  boundary  formulation. 
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Table  1.  Observed  and  model-predicted  amplitudes  (cm.)  and  phases  (degrees,  relative  to  GMT)  for 
the  tide  at  Galveston,  Texas;  Port  Aransas,  Texas;  and  South  West  Pass,  Louisiana.  The  model 
forcing  used  the  standard  Reid  and  Bodine  boundary  formulation. 


Model  Observed 


Station 

Amplitude 

Phase 

Amplitude 

Phase 

Galveston 

9.1 

259.8 

13.5 

275.1 

Port  Aransas 

6.2 

238.0 

7.7 

262.3 

South  West  Pasi 

0.2 

333.5 

1.7 

127.2 

Table  2.  Observed  and  model-predicted  amplitudes  (cm.)  and  phases  (degrees,  relative  to  GMT)  for 
the  M2  tide  at  Galveston,  Texas;  Port  Aransas,  Texas;  and  South  West  Pass,  Louisiana.  The  model 
forcing  used  the  optimized  Reid  and  Bodine  boundary  formulation. 


Model  Observed 


Station 

Amplitude 

Phase 

Amplitude 

Phase 

Galveston 

13.2 

262.1 

13.5 

275.1 

Port  Aransas 

7.9 

242.0 

7.7 

262.3 

South  West  Pass 

3.1 

136.1 

1.7 

127.2 

Table  3.  Observed  and  model-predicted  amplitudes  (cm.)  and  phases  (degrees,  relative  to  GMT)  for 
the  M2  tide  at  Galveston,  Texas;  Port  Aransas,  Texas;  and  South  West  Pass,  Louisiana.  The  model 
forcing  used  the  optimized  Reid  and  Bodine  boundary  formulation  for  the  eastern  open  boundary  and 
again  for  the  southern  open  boundary. 


Model  Observed 


Station 

Amplitude 

Phase 

Amplitude 

Phase 

Galveston 

15.4 

251.0 

13.5 

275.1 

Port  Aransas 

9.6 

234.2 

7.7 

262.3 

South  West  Pass 

8.1 

120.5 

1.7 

127.2 

the  errors,  we  calculated  a  weighted  average  per¬ 
cent  error  for  these  six  northern  Adriatic  Sea  sta¬ 
tions.  The  weighting  is  the  respective  tidal  am¬ 
plitudes  at  each  station.  Thus,  the  errors  at 
stations  with  relatively  small  tidal  amplitudes  do 
not  influence  the  average  error  as  much  as  errors 
at  stations  with  larger  amplitudes.  The  average 
weighted  error  for  the  prediction  of  amplitudes 
using  the  finer  grid  model  is  about  16%.  The 
prediction  of  the  M2  tides  in  the  northern  Adri¬ 
atic  using  the  coarser  grid  model  (Fig.  6)  are 
considerably  worse,  with  percent  errors  reaching 
as  high  as  166%.  It  is  clear  that  the  M2  tidal  pre¬ 
diction  in  the  northern  Adriatic  using  this  coarser 
resolution  model  is  relatively  poor. 


To  simulate  the  tides  in  the  northern  part  of 
the  Adriatic  Sea  only,  a  third,  limited  area  model 
(LAM)  was  produced.  The  LAM  was  based  on 
a  portion  of  the  mesh  of  the  finer  grid  Mediter¬ 
ranean  Sea  model  (same  grid  location  and  reso¬ 
lution;  see  Fig.  7).  The  open  boundary  for  the 
LAM  was  placed  between  points  with  coordinates 
42.53°N,  13.57°E  and  44.56°N,  14.32°E  (Fig.  7). 
The  bathymetry  for  the  LAM  ranged  from  2.3  m 
to  80  m.  The  LAM  was  coupled  with  the  coarse 
grid  model  by  using  scheme  1  and  using  the  op¬ 
timized  Reid  and  Bodine  OBC  in  (13). 

In  the  first  coupling  experiment,  the  ampli¬ 
tudes  and  phases  of  the  sea  surface  elevation  from 
the  coarse  grid  simulation  were  interpolated  to 
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the  open  boundary  of  the  LAM  and  used  as  r/j  in 
13).  We  conducted  coupling  tests  with  Xt  =  1  in 
13)  (standard  Reid  and  Bodme)  and  with  our  op¬ 
timized  approach  where  A*  is  calculated  for  each 
time  step.  The  results  are  shown  in  Fig.  8.  One 
can  see  that  the  use  of  the  standard  Reid  and  Bo- 
dine  boundary  condition  underestimated  the  am¬ 
plitudes  by  up  to  50%.  The  results  of  simulations 
with  the  optimized  OBC  showed  better  agree¬ 
ment  with  the  observations,  but  some  amplitudes 
were  overestimated.  The  average  weighted  error 
was  around  10%,  relatively  good  considering  that 
the  coarse  resolution  model  output  was  used  to 
fr'ce  the  model. 

In  our  experiments,  we  placed  the  open 
boundary  of  the  LAM  in  such  a  way  that  another 
tidal  station,  Porto  Cigale  (44.32°N,  14.27° E), 
lies  on  the  open  boundary.  For  the  second  exper¬ 
iment,  the  amplitude  and  phase  of  t]t  were  equal 
to  the  amplitude  and  phase  of  the  Porto  Cigale 
station  for  each  grid  cell  on  the  open  boundary  of 
the  LAM.  The  results  of  the  simulation  are  shown 
in  Fig.  9  for  the  Reid  and  Bodine  and  optimized 
open  boundary  conditions,  respectively.  The  re¬ 
sults  are  inferior  for  both  open  boundary  condi¬ 
tions  compared  to  the  first  experiment.  Obvi¬ 
ously,  the  observational  data  of  the  one  station  is 
not  enough  to  force  the  limited  area  model  prop¬ 
erly. 

For  the  third  coupling  experiment,  we  com¬ 
bined  the  output  of  the  coarse  grid  run  and  the 
data  of  the  Porto  Cigale  station  to  create  tjt-  For 
the  open  boundary  points  close  to  Porto  Cigale, 
the  values  of  tjt  were  equal  to  that  station’s  data. 
For  the  remaining  grid  points,  the  values  of  tjt 
were  the  same  as  in  first  experiment;  i.e.,  inter¬ 
polated  from  the  output  of  the  coarser  grid  sim¬ 
ulation.  The  results  of  the  LAM  simulations  are 
shown  in  Fig.  10.  Again,  the  standard  Reid  and 
Bodine  boundary  condition  produced  percent  er¬ 
rors  on  the  order  of  50%  in  the  prediction  of  am¬ 
plitudes.  Using  the  optimization  technique,  the 
LAM  predictions  had  a  weighted  average  error  of 
10%  for  the  amplitudes  and  1.7%  for  the  phases. 
It  would  appear  that  forcing  the  model  with  op¬ 
timized  OBCs  and  using  both  observations  and 
the  results  from  the  coarser  model  simulation  pro¬ 
duces  results  that  are  even  better  than  those  of 
the  finer  grid  model  of  the  entire  Mediterranean 
Sea  (Fig.  6).  Additional  research  of  this  phenom¬ 
ena  is  being  pursued. 


7.  CONCLUSIONS  AND  FUTURE  PLANS 

As  we  have  shown,  many  of  the  familiar  for¬ 
mulations  used  as  open  boundary  conditions  in 
barotropic  numerical  simulations  can  have  certain 
optimization  properties  bused  on  the  coefficient 
At.  This  coefficient  allows  the  adaptation  of  the 
boundary  condition  to  the  change  in  the  flux  of 
energy  penetrating  the  open  boundary  as  well  as 
the  minimization  of  differences  between  observa¬ 
tions  and  predictions.  By  choosing  the  appropri¬ 
ate  functional  to  minimize  along  the  open  bound¬ 
ary,  the  proposed  approach  allows  a  modeler  to 
generate  different  types  of  boundary  conditions 
based  on  a  priori  information  and  the  inclina¬ 
tions  of  the  modeler.  Optimized  OBCs  can  be 
considered  as  data  assimilation  schemes,  where 
functions  tjt  ,  ur(s,<)  and  Pt  include  a  prion 
information  of  the  phenomenon  being  modeled. 

The  results  of  simulations  showed  that  the 
proposed  optimized  OBCs  for  a  barotropic  model 
work  very  well  in  reproducing  tidal  phenomenon 
and  can  be  used  in  coupling  coarse  resolution, 
basin  scale  models  with  fine  resolution,  limited 
area  models.  We  plan  to  extend  our  approach  to 
the  modeling  of  three-dimensional  open  bound¬ 
ary  conditions  and  developing  techniques  for  cou¬ 
pling  basin  and  coastal  models  by  using  optimized 
OBCs. 
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Fig.  10.  Same  as  Fig.  8  but  the  model  forcing  used  the  interpolated  output  from  the  coarser  grid 
Mediterranean  Sea  model  plus  the  observed  data  from  Porto  Cigale  station. 
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